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We discuss the formation of quarkyonic chiral spirals in the presence of a magnetic field. The 
explicit breaking of the rotational symmetry by the external magnetic field gives rise to an additional 
chiral spiral that varies along the field direction and rotates in the chiral space between pion and 
magnetic moment components. 

I. INTRODUCTION 

Quarkyonic matter (QyM) is a large N c phase of cold dense quark matter recently suggested by McLerran and 
Pisarski [T] . The main feature of QyM is the existence of asymptotically free quarks deep in the Fermi sea and confined 
excitations at the Fermi surface. The quarks lying deep in the Fermi sea are weakly interacting because they are 
hard to be excited due to Pauli blocking. Their interactions are hence very energetic and the confining part of the 
interaction does not play any role [5] • On the other hand, excitations of quarks within a shell of width Aqcd from 
the Fermi surface interact through infrared singular gluons at large N c and hence are confined. 

In has been recently argued 3 , 4 J that chiral symmetry can be broken in QyM through the formation of a transla- 
tional non-invariant condensate that arises from the pairing between a quark with momentum p, and the hole formed 
by removing a quark with opposite momentum p from the Fermi surface. The inhomogeneous condensate that forms 
in QyM is a linear combination of the chiral condensate (ifnp), and a spin-one, isosinglet odd-parity condensate of 
(■0(T Oz vb). Here z is the direction of motion of the wave. At each given patch of the Fermi surface, z is the direction 
perpendicular to the surface. The authors of Ref. [3] called this combination of two inhomogeneous condensates 
a Quarkyonic Chiral Spiral (QyCS). The (ipu 0z ip) component corresponds to the condensation of an electric dipole 
moment. The QyCS is then characterized by a spatial oscillation between chiral and electric dipole condensates that 
breaks parity and gives rise to an inhomogeneous electric field. 

On the other hand, a common feature of heavy-ion collisions is the generation of strong magnetic fields that are 
produced in peripheral collisions by the positively charged ions moving at almost the speed of light. For the Au-Au 
collisions at the Relativistic Heavy Ion Collider (RHIC) at BNL the field produced is estimated to be ~ 10 19 G [5]. 
Even though this magnetic field decays quickly, it only decays to a tenth of the original value for a time scale of order 
of the inverse of the saturation scale at RHIC |5| , hence it may influence the properties of the QCD phases probed 
by the experiment. Strong magnetic fields will likely be also generated in future experiments planned at the Facility 
for Antiproton and Ion Research (FAIR) at GSI, the Nuclotron-Ion Collider Facility (NICA) at JINR, and the Japan 
Proton Accelerator Research Complex (JPARK) at JAERI, all of which intend to complement the experiments at 
RHIC by reaching regions of even higher densities and intermediate to low temperatures in the QCD phase map. 

In this paper we discuss the effects of an external magnetic field on QyM. As will be seen below, the presence of a 
magnetic field of strength comparable to the square of the QCD scale gives rise to the formation of a second chiral 
spiral, given by a spatial oscillation between a pion condensate {ip^^ip) and a spin-one condensate (^7 7 The 
spin-one condensate corresponds to an inhomogeneous magnetic moment in the direction of the field. 

II. QUARK SELF-ENERGY IN A MAGNETIC FIELD 

We are interested in studying the large N c properties of dense quark matter in the presence of a magnetic field B 
within a region where screening effects are not strong enough to eliminate confining and hence the quarkyonic phase 
can exist. Such a region can be defined by the condition mo <C Aqcd <C fJ., with mo the screening mass and fj, the 
quark chemical potential [1]. 

Our goal is to investigate the self-energy for quarks near the Fermi surface and in the presence of an external 
magnetic field. With that aim, we need to consider the Schwinger-Dyson equation for the quark self-energy Tib(x, x') 
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in a magnetic field, 

Eb(x, x') = g 2 ^t A G B (x, x')j"t B D AB (x - x') (1) 

with Gb(x,x') the full quark propagator in the presence of B and D A ^(x — x') the gluon propagator. Without loss 
of generality we choose the gauge = (0,0, —Bxi,0) for the external electromagnetic field. It produces a constant 
magnetic field in the ^-direction. 

Using Ritus's method [BJ, the transformation to momentum space of the quark propagator and self-energy can be 
done with the help of the matrix eigenf unctions 

E P = ]T E pa (x)A(a), (2) 

tr=±l 

of the asymptotic states of the charged fermions in a uniform magnetic field, where E pcr (x) = 
(47r|e/B|) — e i(p Q x +p 2 x +p 3 x )£) n (p) are the corresponding eigenf unctions, with D n (p) the parabolic cylinder functions 

of argument p = y/2\ejB~\(x 1 — p 2 /efB), and A(er) = /+ '" T 2 7 7 being the spin projectors. The E p<7 (x) depend on the 
spin projection a and the Landau level / > through the non-negative integer n — nil, a) = I + sgn(e/i?)| — |. e/ 
is the electric charge of flavor /. 

In momentum space Eq. ([I]) becomes 



{2^r5\p^p')^{p)Ii l = gh A t B J j^D A u B (q) 



drxe- iq - x ¥, p {xWE p „{x)Ui„G l [p",p) 



A'e'^V^)?"^'^), (3) 



where 



and we used 



and 



E(p,p') = J d 4 xd 4 a;'Ep(a;)S(a;,a;')Ep(a;') 

= (2n) i 5^(p-p')U l j: l (p), (4) 



I-^ U U- ,') -. I 0- )e -^-^D^(q), (5) 



G B (x,x') = ^^^ p „( x )U vl G l (f,p)E p „(x'), (6) 

Here [G'(p)]- 1 = [ 7£ ,P e + ^(IpIT 1 and E'Q?) = E„ (\p\h\\P 11 + ^(Mh^ + E m (|p|), with ^ = 
(poj 0, sgn{efB)^/2\efB\l,pz). The symbol || indicates longitudinal components (0,3), while _L indicates the com- 
ponents perpendicular to the external magnetic field. G l (p,p) is defined by G l (p) with the replacement po — > po + p- 
In the above equations we used the following notation ^ji 4 p = Y^iJ dpadp 2 dpz, E p (x) = 7°E p 7°, and 5^(p — p') = 
8 U 5(po — p' Q )5(p 2 — p' 2 )S{p3 — p' 3 ). The projector II; = 5°' A(sgn(e/i?)) + 1(1 — 5 01 ) appears due to the lack of spin 
degeneracy of the lowest Landau level (LLL) [TJ. 

To consider confined excitations near the Fermi surface, we can use the Gribov-Zwanziger gluon propagator [8], 



D iu (g) = -i5ab ^J f ,^, 2 go»gov, (7) 



where Cf = (N 2 — l)/2N c and the string tension r ~ ^qcd- 

Since the strength of the magnetic fields produced in the heavy-ion collisions is comparable to the QCD scale, it 
makes sense to assume that the external magnetic field satisfies efB < 0(Aq CD ). Summing in the color indexes and 
following derivations similar to those performed in [H| , we can integrate in x, x' and p" , and then transform to polar 
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variables q± = ^fq^ 
indexes to arrive at 



^2, <f = arctan((/2 / Qi) , integrate in ip, and finally perform the sums in Landau levels and spin 



2 8^r 



n, 



d 2 q\\dqj _^ Q G l (p-q, fx)j 
(2tt) 3 



(q 2 ) 2 



(8) 



?l+?2 

2\e f B\- 



where p - q = (p - qo, 0, sgn(e f B)^/2\e]l3\l,p 3 - q 3 ) and q\ 

Because of the Landau quantization of the transverse momentum, the Fermi surface now is formed by a discrete 
set of circles, obtained from the intersection of the cylinders of radius yj2\efB\l with the spherical Fermi surface of 
the system at zero field, as represented in Fig. la. 
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FIG. 1: a) Fermi surface in the presence of a magnetic field; b) Patches 



After transforming to Euclidean space, we can readily integrate in q±, taking into account that the fermion propa- 
gator in Eq. (p| does not depend on this variable, 



s ; (p)n ; 



Nc9 2 2D 



IF 



% 



(27T) 



9,M)74-2 

q-i 



(9) 



Here we introduced a two-dimensional coupling constant g\ D = ^jf- ■ Since pi only enters in the rhs of &9 1 in the quark 
propagator, it can be eliminated with a change of variable. This implies that the self energy is actually independent 

of Pi. 



III. MAGNETIC QUARKYONIC CHIRAL SPIRAL 



We are interesting in exploring the contribution to the self-energy from excitations about the Fermi surface on the 
patches of size ~ ^-qcd *C /i, situated on opposite sites of the Fermi surface, and which lie perpendicular to the 
magnetic field direction (north and south pole patches indicated in Fig. lb) . Quark excitations in these patches satisfy 
P3 ~ A* + <5P3; with Sps < Aqcd "C A*j f° r the momentum component parallel to the field, and p± = y / 2\ejB\l < 
Aqcd for the transverse component, so the number of Landau levels that can contribute to the patch is limited to 

A 2 

L = Irnax = [jjjfjgj], with [...] meaning integer part. These conditions, together with the strong infrared behavior 
in 173 of the gluon propagator, allow us to neglect y/2\ejB\l in \p — q\ in the integrand and similarly drop the term 
^±72\/^\ e fB\l in I]' in both sides of the equation. 

As a consequence, any explicit dependence on the Landau level disappears in Eq. (J9J), except for the projector IF. 
If I 7^ 0, the two spin projections contribute and II; becomes the identity matrix. In this case, Eq. ([9| looks exactly 
like the self-energy equation for two-dimensional QCD in the axial gauge [3l Hj . Eq. (|9| is identical to the one found in 
QyM at zero magnetic field [21 13] , but since here we have one equation for each I contained in the patch, the Landau 
level becomes a flavor index. Hence, there will be L identical equations of the form 

£^3) * J ^74G'(i<?4 + M, 93-^3)74^ l = l,2,...L (10) 
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On the other hand, if I = 0, 11; selects one of the two spin projections only, "up" for positively charged quarks 
(e/ > 0), "down" for negatively charged quarks (e/ < 0). 

Then the original 4-dimensional theory for a quark field with electric charge e t in the presence of a magnetic field 
maps into a (l+l)-dimensional theory with flavor symmetry SU(2L)x U(l) described by the Lagrangian 

Llf f = $o[ir"(d M + ig 2D A, l + r°/i]$o 

L - 1 

+ J2HW(d fl + ig 2 DA lJt )+T f 4$i--trGl„ (11) 
i=i 

Here the spinor fields are defined by <E>q = (ip^\o) and $f = (cp^P ,ip^p), with j" and J, two flavors in the (1+1)-D 
theory that correspond to the up and down spin components of the spinor in the 4D theory. The 2D Dirac T matrices 
are defined in term of the Pauli matrices as T° = a 1 ; T z = —icr 2 ; T 5 = a 3 . 
We can now perform the transformation of the quark fields 

$ ; = exp{-inzT 5 )^[ I = 0, ...L (12) 

to eliminate the chemical potential, which actually remains in the theory through the anomaly of the baryon charge, 
and obtain 

1 



L eff = WC 5 * + ig* D AJ + rV]$; - -TytrG 1 

L 



2"' 

1=0 



J2 <^{d» + iganAJW - \trG\ v (13) 



1=0 



Following arguments along the lines of those discussed in Ref. [3], we can argued that in the above theory chiral 
symmetry is broken through the formation of a chiral condensate for each flavor I, givin g rise to a total chiral condensate 



of the form ($ $') = X^o^i"^;)- Performing the inverse transformation of (12), we obtain two inhomogeneous 



condensates, oscillating in z, and with equal amplitude forming a spiral in the chiraT space: a quarkyonic chiral spiral, 
very similar to the one found at zero field 

= cos(2^z) (IV) ($r 5 $) = -ism(2(jtz) ($V) (14) 
One can easily check that the amplitude of the chiral spiral found in [3] can be expressed as a sum of condensates 
of spin-up and spin-down flavors, ($ $') = (^^f) + (^4.^4.) ■ ^ i s now clear why a condensate of the form r 3 $') = 
l^p'^ipL) — (tp'^tpV) cannot be present in the QyM at zero magnetic field (see Ref. [3] for the definition of the flavor 
matrix r 3 in the present context). There is nothing in the theory at zero magnetic field that can distinguish between 
a condensate of spin-up fields and one of spin-down fields, so physically these two condensates should be the same 
and they must cancel out in ($ 73$'), in agreement with the results of Ref. [3]. 

However, a magnetic field can change the value of ($ t 3 $') because the LLL flavor in the 2D theory has only one 
spin-flavor component. In consequence, two independent chiral condensates can be formed with the primed fields. In 
addition to a chiral condensate of the form 

<$V) = ( ^v o) > + x>; ( v°> + <?iv } >] 
1=1 

we also have now 

1=1 

which is always different from zero due to the LLL contribution, even if the two spin-flavor terms in the sum cancel 
out. The new condensate gives rise to the following second chiral spiral 



($r 3 $) = cos(2^z)($ r 3 $') ($r 3 r 5 $) = -isin(2^z)($ t 3 $') (17) 
: fields in the (3+l)-dimensional theory, the two chiral spirals in the presence of a magnetic 

(fii/j) = Ai cos(2^z) (Vn'V?/') = Ai sin(2^z) (18) 



(V'7 1 7 2 V') = A 2 cos(2^z) (ipj 5 ip) = A 2 sia(2/zz) 



(19) 
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IV. CONCLUSIONS 



In this paper we have considered QyM in a magnetic field of strength close to, but not larger than the QCD scale. 
We showed that in leading order in the large N c limit, a (3+l)-dimensional QCD theory with one flavor can be 
mapped into a (l+l)-dimensional QCD theory with 2L + 1 flavors and global symmetry SU(2L)x U(l). 

Due to the contribution of the LLL, two independent chiral spirals are formed in the Fermi surface's patches 
perpendicular to the magnetic field. One, similar to the spiral condensate that forms at zero field, is a combination 
of a chiral condensate and an electric dipole condensate, rotating in the chiral space and varying along the direction 
parallel to the field. The other chiral spiral is a combination of a condensate of magnetic moment and a pion 
condensate, also varying in the direction parallel to the field. Both parity and time-reversal symmetries are broken 
in the system. 

The spontaneous generation of inhomogeneous condensates with electric and magnetic dipole moments may lead to 
interesting observational implications. A study of those potential consequences in dense environments like the cores 
of neutron stars or the planed high-density heavy-ion collisions experiments will be addressed in future works. 
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